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ABSTRACT 

Here we present timing analysis of a set of eight Rossi X-ray Timing Explorer (RXTE) observations of the 
black hole candidate GX 339-4 that were taken during its hard/low state. On long time scales, the RXTE All Sky 
Monitor data reveal evidence of a 240 day periodicity, comparable to timescales expected from warped, precessing 
accretion disks. On short timescales all observations save one show evidence of a persistent /qpo ~ 0.3 Hz QPO. 
The broad band (10~ 3 -10 2 Hz) power appears to be dominated by two independent processes that can be modeled 
as very broad Lorentzians with 2^1- The coherence function between soft and hard photon variability shows 
that if these are truly independent processes, then they are individually coherent, but they are incoherent with one 
another. This is evidenced by the fact that the coherence function between the hard and soft variability is near 
unity between 5 x 10~ 3 -10Hz but shows evidence of a dip at / m 1 Hz. This is the region of overlap between 
the broad Lorentzian fits to the PSD. Similar to Cyg X-l, the coherence also drops dramatically at frequencies 
> 10Hz. Also similar to Cyg X-l, the hard photon variability is seen to lag the soft photon variability with the 
lag time increasing with decreasing Fourier frequency. The magnitude of this time lag appears to be positively 
correlated with the flux of GX 339-4. We discuss all of these observations in light of current theoretical models 
of both black hole spectra and temporal variability. 

Subject headings: accretion — black hole physics — Stars: binaries — X-rays: Stars 



1. INTRODUCTION 



analysis in general, has been presented by van der Klis (1989) 



In a companion paper to this work (Wilms et al. 1998a; here- 
after paper I) we have presented spectral analysis of a series of 
Advanced Satellite for Cosmology and Astrophysics (ASCA) 
and simultaneous radio/Rossi X-ray Timing Explorer (RXTE) 
observations of the black hole candidate (BHC) GX 339- 
4. This source has e x hibited both spectrally soft states (cf. 
Grebenev et al. 1991; Miyamoto et al. 1991) and spectrally 
hard states (cf. [ Grebenev et al. 1991 ; Miyamoto et al. 1992 



Here we apply these Fourier analysis techniques in the same 
manner as for our RXTE observations of Cyg X-l (N owak 
et al. 1998a). Specifically, we u sed the same techn iques for 
estimating: deadtime corrections dZhang et al. 1995fc Z hang & 
Jahoda 1996); the error bars a nd Poisson noise levels of the 
Pow er Spectral Density (PSD) ([Leahy et al. 1983j v an der Klis 



1989); the error bars and noise levels for the coherence function 



(Bendat&Piersol 1986 



Zdziarski et al. 1998). Both the ASCA and RXTE observations 



ror bars and noise leve 



Vaughan & Nowak 1997); and the er 



presented in paper I showed GX 339-4 to be in spectrally hard 
and low luminosity (3-9 keV flux < 10~ 9 ergs cm" 2 s" 1 ) states. 
The eight RXTE observations spanned roughly a factor of five 
in terms of observed 3-9 keV flux. 

In this work we shall consider the timing analysis of the 
eight RXTE observations. Timing analysis for GX 339^4 pre- 
viously has been prese nted for the soft, bright 'very high state' 
( Miyamoto et al. 199 1|), the soft, fainter 'high state' (Grebenev 
et al. 1991), and the hard 'low state' ( prebenev et al. 1991 ' 



s for the Fourier frequency-dependent 
time lag bet ween hard and soft photon variability (B endat & 
Piersol 1986; |Nowak et al. 1998a[ ). A self-contained discussion 
of these techniques as re gards RXTE timing analysis is given 



by Nowak et al. (1998a) 



Miyamoto et al. 1992| ). In addition, timing analysis has been 



This paper is organized as follows. First we consider evi- 
dence for long term periodicities in the RXTE All Sky Monitor 
(ASM) data in section |^. In section |5] we consider the Power 
Spectral Density (PS D). In section ^ we consi der both the co- 
herence function (cf . |Vaughan & Nowak 1997 ) and the Fourier 
frequency-dependent time lags between hard and soft photon 
variability. We discuss the theoretical implications of these ob- 
servations in section ||. We then summarize our results in sec- 
tion^ 

2. ALL SKY MONITOR OBSERVATIONS OF LONG 
'CHARACTERISTIC TIMESCALES' IN GX 339-4 

We used data from the All Sky Monitor on RXTE to study 
the long-term behavior of GX 339-4. The ASM is an array of 
three shadow cameras combined with position sensitive propor- 
tional counters that provides for a quasi-continuous coverage of 
the sky visible from RXTE. In practice, lightcurves in three en- 
ergy bands — 1.3-3.0 keV, 3.0-5.0keV, and 5.0-12.2 keV — as 
well as over the whole ASM band are publically available from 

'JILA, University of Colorado, Campus Box 440, Boulder, CO 80309-0440, USA; {mnowak, dove) @rocinante. colorado.edu 
2 Institut filr Astronomie und Astrophysik, Abt. Astronomie, Waldhauser Str. 64, D-72076 Tubingen, Germany; wilms@astro.uni-tuebingen.de 
3 also, CASA, University of Colorado, Campus Box 389, Boulder, CO 80309-0389, USA 

1 



presented for so-called 'intermediate states', which are bright 
but spectrally harder than the t ypical 'very high' or 'high states' 
(Mendez & van der Klis 1997 ). In general, the harder states ex- 
hibit more te mporal variability than the softer states (cf. van 
der Klis 1989). Timing analyses of the GX 339-4 hard state 
have shown similar resul ts to analyses of other hard state BHC 
( [Miyamoto et al. 1992[ ). Specifically, the GX 339-4 tem- 
poral variabil ity has been observed to b e similar to that of 
Cygn us X-l ([Belloni & Hasinger 19904 Belloni & Hasinger 
1990b; Miyamoto et al. 1992| ; |Nowak et al. 19984 and refer- 
ences therein). 

A discussion of Fourier techniques in specific, and timing 
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Fig. I.— Left: RXTE All Sky Monitor data for GX 339-4 (5 day averages in the 1.3-12.2 keV band) vs. Truncated Julian Date (TJD) = Julian 
Date (JD) -2450000.5. Dashes indicate dates of our RXTE pointed observations, diamonds indicate dates of MOST observations (Hannikainnen 
et al. 1998), and triangles indicate dates of ATCA radio observations (Corbel et al. 1998). Right: A Lomb-Scargle Periodogram (cf., Lomb 1976, 
Scargle 1982) of the ASM data for TJD < 800. We have used 600 periods ranging from 2 weeks to 500 days. Lines are estimates of the 99.9% 
and 90% significance levels. 



the ASM data archives ( Lochner & Remillard 1997 ). Typically 
there are several 90 s measurements available for each day. Fur- 
ther descr iptions of the instrum ent and first year results are pre - 
sented by Levine et al. (1996)| and Remillard & Levine (1997) . 

In Figure [l] we present the ASM data of GX 339-4 up un- 
til Truncated Julian Date (TJD) « 800 (1997 November 15). 
See paper I for a presentation of the complete to date ASM 
data of GX 339-4. We also indicate in this figure the dates 
of our RXTE observations, as well as the dates of Australian 
Telescope Compact Array (ATCA) and Molongolo Observing 
Synthesis Telescope (MOST) radio obse rvations of GX339-4 
( Render et al. 1997|; £ orbel et al. 1997t H annikainnen et al 



parent in the ASM data during the hard stat e. A 198 day period 



icity also has been o bserved in LMC X-3 ( Cowley et al. 1991 
Wilms et al. 1998q ). We note that the observed timescales of 



1998; paper I). The average la error bar for a point presented 
in Figure]l]is w 0.4 cts s" 1 . Note the three peaks (at TJD w 280, 
480-580, 780) in the ASM data that occur at intervals separated 
by w 250 days. 

We determined the significance of any possible long term pe- 
riodicities in the AS M light curves by computin g the Lomb- 
Scargle Periodogram dLomb 1976| ; |Scargle 1982] ) for the 1.3- 
12.2 keV band for 5-day average data. Here we only average 
data where the best fit to the sourc e position and flux in an ASM 
observation has a x%d < 1-5 (cf. ^ochner & Remillard 199% ) 
in each of the three ASM energy channels. We only considered 
ASM data taken before TJD 800, as GX 339-4 underwent a 
state change shortly thereafter (see paper I). The periodogram 
shown in Figure [y reveals evidence of a 240 day period at 
greater than the 99.9% significance level. [The significance lev- 
els we re estimated following the methods of H ome & Baliunas 
(1986)TEp och folding (cf ILeahy et al. 19831; Schwarzenberg- 
Czerny 1989; pavies 1990| ) of the ASM lightcurves also shows 
evidence of this 240 day periodicity. 

Long timescale periodicities and quasi-periodicities are rela- 
tively common in ASM observations of binary sources (Remil- 
lard 1997, Private Communication). Evidence for a 294 d pe- 
riodicity i n Cygnus X-l has been previously reported ( Kemp 
et al. 1983; |Priedhorsky, Terrell & Holt 1983| ), and is readily ap- 



these periodicities are comparable to the timescales expected 
from the the radiation press ure driven warping instability dis- 
covered by |Pringle (1996)1 (see also Ma loney, Begelman & 
Pringle 1996; Maloney & Begelman 19971; M aloney, Begelman 
& Nowak 1998). This is a fairly generic instability that causes 
a radiatively efficient (i.e., non-advection dominated) accretion 
disk to warp and precess on O(100 day) timescales. Such a 
timescale is consistent with the periodicities seen in many bi- 
naries and with the 240 periodicity that we see in GX 339—4. 

As discussed in paper I, however, both the Compton corona 
models and the ADAF models of the spectral data suggest that 
the observed flux variations are in large part attributable to vari- 
ations of the coronal radius. A pure warped, precessing disk 
model would invoke only inclination angle effects. This is un- 
likely to be the case for G X 339-4. As discussed by M aloney, 
Begelman & Pringle (1996), however, the precession timescale 
and warp shape are more sensitive to the outer boundary condi- 
tions than the inner boundary conditions. Perhaps then it is pos- 
sible that what we observe is a combination of a quasi-steadily 
precessing disk on large radii combined with coronal structure 
changes on small radii. Although 240 days is a characteris- 
tic timescale in both the Lomb-Scargle periodogram and in the 
epoch folding analysis, it is obvious from Figure [j] that we are 
not observing a strictly periodic phenomenon. 

3. POWER SPECTRAL DENSITIES 

We studied the variability of GX 339—4 using data from 
the Proportional Counter Array (PCA) on-board RXTE. The 
PCA consists of five nearly identical co-aligned Xenon propor- 
tional counter units (PCUs) with a total effective area of about 
6500 cm 2 , and it is sensitive in the energy range from 2 keV to 
~ 60keV (lahoda et al. 1996). We only used data where all five 
PCUs were turned on, and where the elevation angle between 
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the spacecraft pointing direction and the limb of the earth was 
greater than 10°. PCA count rates for the GX 339-4 obser- 
vations ranged from 200 to 800 cts s , and were of O(10ks) 
duration. 

To study the short timescale variability, we created 2" 8 s res- 
olution lightcurves in 4 different energy bands: 0-3.9 keV, 3.9- 
7.5keV, 7.5-10.8keV, 10.8-21.9 keV; hereafter labelled bands 
A, B, C, and D 1 . Note that this is one fewer energy band and 
half the time resolu tion than for our Cyg X-l observations 
(Nowak et al. 1998a). This was required to obtain good statis- 
tics because the GX 339—4 observations discussed here ranged 
from 200-800 cts s , as opposed to the 4500 cts s" 1 observed 
for Cyg X-l. Energy bands A-D each had roughly the same 
count rate for a given observation. As for the spectral analy- 
sis (paper I), we found that all the observations were similar in 
their properties, with the exception of Observation 5. Observa- 
tion 5 was approximately a factor of 5 fainter than the brightest 
observation (Observation 1). 

As for our RXTE observations of Cyg X-l ( Nowak et al. 
1998a), we combined low frequency («8x 10~*-7 x 10" 3 Hz) 
FFTs made from a few (6-16) data segments of length 512— 
1024 s with higher frequency (7 x 10~ 3 -128Hz) FFTs made 
from somewhat more (^ 70) data segments of length 128 s. 
From these data segments we created PSDs for all of our ob- 
servations. We present two examples of the resulting PSDs 
in Figure ||. With the exception of Observation 5, all PSDs 
were qualitatively and quantitatively similar to that presented 
in Figure H for Observation 7. The PSDs for Observations 1-4 
and 6-7 have shapes roughly similar to what we observed for 
Cyg X-l: flat from « 10~ 3 -0.03Hz, approximately oc in- 
betwe en 0.03-3 Hz, and steeper above « 3 Hz ( Nowak et al. 
1998a). The PSDs for Observation 5 have qualitatively simi- 
lar shapes; however, the break frequencies are approximately a 
factor of three lower. 

Doubly broken power laws provided adequate descriptions of 
the P SDs for our RXTE observations of Cyg X-l (N owak et al. 
1998a). Here, however, we clearly see that the GX 339-4 PSDs 
contain more structure. An w 0.3 Hz quasi-periodic oscillation 
(QPO) is evident in Observation 7. In fact, this QPO appears 
at some level in all of our observations except Observation 5. 
Even ignoring the QPO features, the PSD is only marginally 
approximated by a doubly broken power law. 

We have tried to fit the shape of the PSDs with a number of 
different functional forms. For example, a singly broken power 
law with a Lorentzian absorption feature at « 0.3-1 Hz pro- 
vides good fits to all of the PSDs. Here we show the results of 
fitting a weak power law plus two broad Lorentzians. As we 
further elaborate upon in §Q and §|J such a fit may have some 
physical relevance. The Lorentzian and power law components 
may represent separate, broad-band processes that are individ- 
ually coherent (between their soft and hard variability), but that 
are incoherent with one another. 

Sample fits to the broad-band power are presented in Fig- 
ure U and results for all of our data are presented in Table [j]. 
The functional form of the Lorentzians that we fit is given by 



R 2 Qh 



ti + Q 2 (f-fo) 2 



(i) 



maximum of the Lorentzian), and R is the fit amplitude (root 
mean square variability, rms = fl[l/2-tan- 1 (-0/7rj I/2 , i.e. rms 
= R as Q — > oo). For all the fits presented in Table |J Q ^ 1 . We 
thus cannot consider the broad-peaks in the observed PSDs to 
be "quasi-periodic oscillations". The fits are more rightly con- 
sidered to be indicative of broad-band power with characteristic 
timescales ~ /y 1 . Note also that, for all of the fits, x 2 e d ~ 2- 
3. Although this formally represents an unacceptable fit, this 
is partly attributable to the excellent statistics achievable with 
RXTE. Even small fractional deviations from the fit are highly 
statistically significant. In practice we have found that x 2 ed ^ 2 
is nearly impossible to achieve with any simple functional fits. 
Futhermore, the fits with the largest x 2 e d showed the most pro- 
nounced 0.3 Hz QPO, as in Observation 7. 

For Observations 1-4 and 6-7, the low frequency broad-band 
power was peaked near « 0.3 Hz, while the high frequency 
broad-band power was peaked near ss 2.5 Hz. The rms vari- 
ability amplitudes were all comparable for these observations. 
Observation 5 had a somewhat larger rms variability amplitude. 
The trend for all observations was for the rms variability ampli- 
tude, measured from 7 x 10~ 3 ^40Hz, to decrease with increas- 
ing energy band, although this is only marginally evident in the 
individual fit components (cf. Figure [j]). In Figure [j] we plot 
the best fit rms variabilities and peak frequencies for the two 
Lorentzian fit components vs. the observed 3-9 keV flux. 

Observations 1-4 and 6-7 span a factor of two in observed 
3-9 keV flux. The PSD parameters, however, show no obvious 
trends with flux, except for exhibiting weak evidence for the 
rms variability of the low frequency Lorentzian fit component 
to increase with 3-9 keV luminosity. Perhaps the most remark- 
able aspect of this subset of the observations is how similar all 
the PSDs appear to one another despite the factor of two spread 
in the observed 3-9 keV flux. 

The w 0.3 Hz QPO is most prominent in Observation 7. Its 
strength and width are somewhat difficult to characterize, how- 
ever, as it is difficult to determine the 'continuum' level of the 
PSD to which it should be compared. The rise of the QPO 
appears sharper than its decline, and typically there is a PSD 
'dip' at « 0.05 Hz before the peak of the QPO. We have tried a 
variety of narrow Lorentzian plus broad-band power fits to the 
features; however, the best fit parameters are highly dependent 
upon the assumed form of the underlying broad-band PSD. Fur- 
thermore, narrow Lorentzian fits do not capture the asymmetric 
shape of the QPO. 

Being unable to find a satisfactory functional fit to the QPO, 
we have chosen to characterize it in the following manner. We 
have measured the location of the PSD 'dip' that occurs be- 
fore the QPO peak, as well as the location of the QPO peak. 
(Factoring in noise fluctuations, the QPO is sharp enough from 
dip to peak that the location of each is likely accurate to better 
than two unaveraged frequency bins, i.e. 0.015 Hz.) We have 
also measured the dip and peak PSD amplitudes. As for the 
fits to the broad-band PSD discussed above, there are no obvi- 
ous trends between QPO parameters and the observed 3-9 keV 
flux. Furthermore, there is no obvious trend between QPO peak 
power and frequency, as shown in Figure & Although the QPO 
frequency is characteristically near 0.3 Hz, it is not at a steady 
frequency and it varies over a range of 0.08 Hz. Peak QPO am- 
plitudes and widths, however, are somewhat less variable. Note 
that the total rms variability in the region of the QPO typically 
is « 5%, therefore the QPO amplitude is at most a few percent. 



Here /o is the resonant frequency of the Lorentzian, Q is 
the quality factor (w /o/A/, where A/ is the full-width-half- 

'We also created 2~ 13 s resolution lightcurves to search for signatures of high frequency features. No evidence for such features was found. Note also that the 
effective lower limit for energy band A is more 2 keV. 
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Fig. 2. — PSDs with associated uncertainties for Observation 5 (left) and for Observation 7 (right). All PSD are for the one-sided normalization 
of Belloni & Hasinger (1990b), where integrating over positive frequencies yields the total mean square variability. Crosses correspond to energy 
band A, squares to energy band B shifted upwards by a factor of 10, triangles to energy band C shifted upward by a factor of 100, and diamonds 
to energy band D shifted upward by a factor of 1000. (See text for the energy range of the bands.) Root mean square (rms) variability values were 
calculated between / = 7 x 10~ 3 -40 Hz. 




Fig. 3. — Left: Diamonds are the PSD for Observation 5 (energy bands A-D summed), normalized as in Belloni & Hasinger (1990b). Right: 
Diamonds are the PSD for Observation 7 (energy bands A-D summed). For both figures, the solid line is the best fit power law plus two broad 
Lorentzians. (Fit range is 0.03-10 Hz for Observation 5, and 0. 1-30 Hz for Observaion 7.) Dashed lines show the individual components of the 
fits. 
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Table 1 

Fits to the PSD of the form: A/^+tT 1 (tfQifi/lff + Qiif-fift+RlQifi/lfl + Qlif - fif}) ■ Energy ranges of bands A-D are 
described in the text. t represents energy bands a-d summed together. all fits were in the range / = 0.1-30 hz (47 
degrees of freedom), except for observation 5 which was fit in the range / = 0.03-10 hz (45 degrees of freedom). 
Errors are the nominal 90% confidence level for one interesting parameter (Ax 2 = 2.71). 
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Fig. 4. — Upper Left: rms variability of the low-frequency broad Lorentzian fit component vs. the 3-9 keV flux of the observation. Upper Right: 
Peak frequency of the low-frequency broad Lorentzian fit component vs. the 3-9 keV flux of the observation. Lower Left: rms variability of the 
high-frequency broad Lorentzian fit component vs. the 3-9 keV flux of the observation. Lower Right: Peak frequency of the high-frequency broad 
Lorentzian fit component vs. the 3-9 keV flux of the observation. Symbols are the same as for Figure 2. 
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4. COHERENCE AND TIME LAGS 



As discussed by Vaughan & Nowak (1997' , the coherence 
function is a measure of the degree of linear correlation between 
two time series. Specifically, it gives the fraction of the mean- 
squared variability at a given Fourier frequency in one time se- 
ries that can be attributed to, or equivalently predicted from, the 
other. The fact that coh erence is often near unity over a wide 
range of frequencies (cf. |Vaughan & Nowak 1997|; N owak et al. 
1998a) is what makes it meaningful to then talk about charac- 
teristic Fourier frequency-dependent time del ays between two 
time series (cf. [Vliyamoto & Kitamoto 1989|; M iyamoto et al. 
1992; |van der Kris 1989| ; ^lowak et al. 1998a| ). 

The techniques that we used to calculate the coherence and 
time delays for the GX 3 39^ 4 lightcurves are disc ussed in 
Vaughan & Nowak (1997)| and |Nowak et al. (1998a)| . For aU 
observations we have calculated the coherence function, 



7 2 (/) : 



| (S* (/)//(/)) I 2 

(IWXIW)! 2 



(2) 



between the FFTs for energy band A [S(f)] and the FFTs for the 
other three energy bands [//(/), angle brackets indicating aver- 
ages over Fourier frequencies and individual da ta segments; cf. 
Vaughan & Nowak 1997|; |Nowak et al. 1998a|], For all obser- 



vations except for Observation 5, we have calculated the time- 
delay between energy band A and the other three energy-bands. 
For Observation 5, noise limitations only allow us to calculate 
the time-delay between the lightcurves for energy band D and 
energy band A. We have averaged over logarithmically spaced 
frequency bins, ranging over frequencies / — > 1 .4/ for all ob- 
servations. 

As shown in Figure ||, the coherence function is near unity 
from 10" 3 to « 3 Hz (Observation 5) or w 10 Hz (all other ob- 
servations). Above w 3-10 Hz, there is a noticeable drop in co- 
herence, simi lar to what we saw with our RXTE observations 
of Cyg X-l ( [Nowak et al. 1998a| ). The Cyg X-l observations 



also showed a loss of coherence below w 0.02 Hz; however, 
the coherence in GX 339-4 remains near unity down to Fourier 
frequencies as low as w 10~ 3 Hz. 

We note that the coherence between energy bands A and D 
shows evidence of an « 0. 1 dip near 0.5 Hz (Observation 5) and 
near 2 Hz (all other observations). The evidence for this dip is 
weak in Observation 5 (comparing band A to D, the 0.45 Hz 
point is 2fj below unity, and the adjacent 0.64 Hz point is only 
la below unity), but is somewhat stronger for Observation 7 
and the other observations (comparing band A to D for Obser- 
vation 7, the 1.8 Hz point is > 3(7 below unity, and the adjacent 
2.5 Hz point is « 2.7ct below unity). These dip frequencies, 
however, are approximately the frequencies at which the two 
broad Lorentzian and power law fit components (cf. Figure || 
and Tab. [j]) overlap. Thus we hypothesize that there are indeed 
multiple broad-band processes occuring in GX 339^1 that are 
individually coherent but that are incoherent with one another. 
As with Cygnus X-l, we further hypothesize that the loss of 
coherence at high-frequency is indicative of mulitple incoher- 
ent process es, possibly 'flares' feeding the corona on dynamical 



timescales ( Nowak etal. 19984 [Nowak et al. 1998bh . 

In Figure ^| we also show the energy-dependent and Fourier 
frequency-dependent time delays. The overall structure is very 



& Kitamoto 1989; 


Miyamoto et al. 1992|; prary et al. 1998 


Nowak et al. 1998; 


). In the regions not dominated by noise 



ton variability, and the delays decrease with increasing Fourier 
frequency. The detailed frequency-dependent structure of the 
delays, however, is more complicated than a simple power law. 
For example, the time lags between bands A and D for Observa- 
tion 7 show a flattened region near 1 Hz. As for Cyg X-l, there 
is a large dynamic range in the time delays, with the longest 
time delays (f» 0.1 s) being much longe r than typical charac - 
teristic timescales of a small corona (cf. Nowak et al. 1998b). 



Also as ha s been observed for Cyg X-l (Miyamoto & Kita- 
moto 1989; [Vliyamoto et al. 1992|; |Crary et al. 1998j; N owak 
et al. 1998a), the time delay observed in GX 339-4 is consis- 
tent with scaling as the logarithm of the ratio of the two energies 
being compared. T his latter fact has prompted K azanas, Hua 
& Titarchuk (1997) (hereafter KHT) to suggest that the time- 
delays are related to photon propagation timescales in a very 
large (R « 10 5 GM/c 2 ) corona. In Figure [7] we show the mea- 
sured time delay for three Fourier frequencies (0. 1 Hz, 0.9 Hz, 
2.5 Hz) as a function of the measured 3-9 keV flux. The low- 
est flux observation consistently shows the shortest time delays 
at nearly all Fourier frequencies. The highest flux observa- 
tion shows the longest time delays at many Fourier frequen- 
cies, and shows at least the second longest time delays at nearly 
all Fourier frequencies. Observations at intermediate fluxes are 
scattered both positively and negatively about an extrapolation 
between the low- and high-flux point. It has been previously 
shown that the time delay de creases in Cyg X-l as it transits 
from the hard to soft state (cf. |Cui et al. 1997 ). Here we present 
a possible correlation between the magnitude of the time delay 
and the energy flux for the hard state within a single source. 
The strong decrease of the time delay for the lowest flux ob- 
servation is counter to the simplest expectations if the coro- 
nal size increase s with decreasing flux (e.g. Esin, McClintock 
& Narayan 1997), or if characteristic 'propagation speeds' in 
the co rona decrease with decreasing luminosity (N owak et al. 
1998b). However, if the coronal size is decreasing with decreas- 
ing luminosity then this observed decrease is understandable in 
terms of propagation models, whether it be propagation of pho- 
t ons (KHT) or propa gation of some sort of other disturbance 
(|Nowaket al. 1998bb. 



statistics, the hard photon variability always lags the soft pho- 



We have found one potential correlation between the mea- 
sured time lags and the coherence function, related to the vec- 
tor analog_y_fOTthejjhase lags and coherence function discussed 
by Nowak et al. (1998a). If we take the Fourier transform of a 
soft X-ray lightcurve, A s (f), and a Fourier transform for the 
hard X-ray lightcurve, A/,(/), the cross spectrum is given by 
A*(/)A;,(/), which can be considere d as a vector in the complex 
plane (cf. Fig. fy. As discussed by Vaughan & Nowak (1997). 
the magnitude and angle (corresponding to the phase delay be- 
tween hard and soft variability) of this vector is well-defined if 
the coherence function is unity. 

If the net observed cross spectrum, however, is made up of 
the sum of individually coherent processes, it is possible that 
the net observed coherence will be less than unity. As was noted 
by Vaughan & Nowak (1997), specifically eq. (10), a sum of 
individually coherent processes is itself coherent if and only if 
each process has the same identical transfer function from soft 
to hard photon variability. In terms of the vector analogy, this 
is a statement that the vectors representing each process within 
the sum all lie in the same direction. The coherence function in 
such a case is the square of the ratio of the magnitude of this 
sum compared to the magnitude of the sum if all the vectors 
had the same angle. 

Let us consider the following special case of eq. (10) from 
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Fig. 6. — Coherence function and time lags of various energy channels [squares: channel B, triangles: channel C (coherence offset by 0.5), 
diamonds: channel D (coherence offset by 1.0)] as compared to energy channel A. Top Left: Coherence for Observation 5. Solid lines correspond 
to unity coherence. Top Right: Time lags for Observation 5. Crosses are where the soft variability lags the hard variability. Solid line is the 
Poisson noise level for channel D time lags. Bottom Left: Coherence for Observation 7. Solid lines correspond to unity coherence. Bottom Right: 
Time lags for Observation 7. Crosses are where the soft variability lags the hard variability. Solid line is the Poisson noise level for channel D 
time lags. 
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Fig. 7. — Time lags at Fourier frequencies 0. 1 Hz (diamonds), 0.9 Hz (triangles), and 2.5 Hz (squares) as a function of measured 3-9 keV flux. 



Vaughan & Nowak (1997). Assume that there are multiple 
input (soft) processes, indexed by i, with Fourier amplitudes 
A' s (f). Let us assume that each of these input processes leads 
to an output process with Fourier amplitude A' h (f) that has a 
coherent phase delay of A.<pt(f). (Here we take the Fourier am- 
plitudes to be real quantities.) Finally, assume that the input 
processes are incoherent with one another, and likewise that the 
output proc esses are incoherent with on e another. Generalizing 
eq. (10) of Vaughan & Nowak (1997), the resulting measured 
coherence function will then be 



E A i 2 E A * : 



(A^-A0,) 2 
2 



(3) 



where we have adopted a small angle approximation. Under the 
same assumptions, the net measured phase lag will be given by 



A0 m (/) i 



E^A Aft 



(4) 



That is, the measured phase lag is simply the weighted mean of 
the individual phase lags. We illustrate this situation in Fig. [| 
For such a model as this, the net observed phase delay is re- 
lated to the net observed coherence function, and both in turn 
are related to the observed soft and hard PSDs. Fit parameters 
are the amplitudes {A' s , A' h ) of the individual components of the 
soft and hard PSDs, and the phase lags (A</>,) between the soft 
and hard variability for each of these components. 

We have searched for such a connection between the Fourier 
phase lag and the coherence function by applying equations (^J) 
and (Q) to the data. We have assumed that there are three com- 
ponents to the PSD: a low and high frequency Lorentzian and a 
power law, as for the fits presented in Table [|. In the fitting pro- 
cedure we assumed that all three A0,(/) were independent of 
Fourier frequency. We refit the PSD data simultaneously while 



fitting the phase lag and coherence data. We searched for the 
minimum of the sum of the \ 2 from the four data sets being 
fit: soft X-ray PSD, hard X-ray PSD, phase lags, and coher- 
ence. We fit the PSDs over the same frequency range and with 
the same frequency binning as in Table [j]; however, we only 
fit the coherence and phase lags over the range / = 0.2^1-Hz. 
This was the frequency range over which the phase lags were 
least affected by Poisson noise, and it avoided the strong loss 
of coherence at high frequency. These strong high frequency 
coherence losses, as we discuss further below, may be due to 
nonlinear p rocesses, rather than be du e to the sum of linear pro- 
cesses (cf. Vaughan & Nowak 1997[ ). The results for these fits 
are presented in Figure ^| and in Table |2| 

Assuming that the PSD, phase lags, and coherence are the re- 
sult of summing three independent processes (a power-law plus 
two Lorentzians), we see that equations ([}]) and (Q) seem to pro- 
vide a rough description of the time lag and coherence data. The 
fits make plausible that there is indeed a deeper underlying con- 
nection between the time lag and coherence data. One feature 
of these fits is notable. Namely, in order to produce coherence 
drops as large as are seen, one needs to add linear processes 
with greatly varying intrinsic time lags. In fact, one process, 
the power law, has nearly no time delay between soft and hard 
variability, whereas the high frequency Lorentzian process is 
seen to require even intrinsically longer time lags than the al- 
ready very long time lags that are measured. 

5. DISCUSSION 

Let us consider these results in light of two models: the 
Comptoniza t ion model of KHT and 'shot noise' models (cf., 
Ten-ell 1972t Sutherland, Weisskopf & Kahn 1978fc Priedhorsky 
et al. 1979; Miyamoto & Kitamoto 1989fc L ochner, Swank & 
Szy mkowiak 1991; |Nowak 1994 hereafter N94; B elloni et al. 
1997; |Poutanen & Fabian 1998| , hereafter PF; and references 
therein). In the former model, a source of white noise is as- 
sumed to be present at the center of the system, which is then 
"filtered" by passage through a Compton corona. The PSD is 
attenuated on those timescales that are shorter than the char- 
acteristic diffusion timescales through the corona. Therefore, 
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Fig. 8. — Vector analogy for phase lags and coherence. Cross power spectra (Aj A\, A\ A\, A\ A\) can be considered as vectors in the complex 
plane. The observed phase lag, A<j) m , is the angle of the vector sum of the individual components of the cross power spectrum. The coherence 
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Fig. 9. — Left: Phase lags between energy band A and D as a function of Fourier frequency for Observation 5 (top) and Observation 7 (bottom). 
(All points are hard variability lagging soft variability.) Right: [1 -7„(/)], where 7^(/) is the measured, noise subtracted coherence function 
between energy band A and D for Observation 5 (top) and Observation 7 (bottom). All data has been logarithmically binned over frequencies 
/ — > 1.4 /. Solid lines are the best fit results for fitting equations of the form of eq. (2) and eq. (3) (see text and Table 2 for further explanation). 
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Table 2 

Best fit Fourier phase lags between energy channels A and D, where we have used equations of the form of eq. (2) 

AND EQ. (3) AND WHERE WE HAVE ASSUMED THREE SEPARATE COMPONENTS TO THE SOFT AND HARD X-RAY PSD. PSD COMPONENTS 

are as in Table 1: a low frequency Lorentzian (11), a power law (pi), and a high frequency Lorentzian (12). PSD data, 

logarithmically binned over frequencies / — > 1.1/, were included in the fitting process with the resulting fit 
parameters always being within the error bars shown in table 1 . phases — a0/i , a(f> p i, a0; 2 — are the fourier phase 
lags between hard and soft variability for each psd fit component individually, and positive phase corresponds to 
hard variability lagging behind soft variability. the phase lags were assumed to be constant as a function of 

Fourier frequency. 



Obs. 


A</> n 


A0„, 


A(j} l2 


X 2 /DoF 


2 

Xred 




(rad) 


(rad) 


(rad) 






1 


0.13 


0.06 


0.49 


301.5/119 


2.53 


2 


0.18 


0.05 


0.26 


283.4/119 


2.38 


3 


0.13 


0.04 


0.28 


269.9/119 


2.27 


4 


0.15 


0.06 


0.66 


195.8/119 


1.65 


5 


0.06 


0.04 


0.25 


132.4/113 


1.17 


6 


0.16 


0.05 


0.38 


217.4/119 


1.82 


7 


0.11 


0.04 


0.27 


309.9/119 


2.60 


8 


0.16 


0.05 


0.41 


198.2/119 


1.67 



the PSD is more attenuated at higher frequencies than at lower 
frequencies, which leads to the characteristic power law PSDs 
that are observed in hard state BHC. Time delays are created 
by the difference in diffusion times through the corona for hard 
and soft photons. Photons that scatter over large radii will have 
their intrinsic high frequency variability wiped out; therefore, 
any observed high frequency variability must be due to photons 
that scattered over short path lengths. High frequency variabil- 
ity thus exhibits short time delays between hard and soft pho- 
ton variability. Low frequency variability potentially can be ob- 
served from photons that have scattered over large path lengths, 
and thus it can exhibit longer time delays between hard and soft 
variability. Extremely large coronal sizes of O(10 5 GM/c 2 ) are 
required to produce the longest observed time lags. 

In shot noise models the lightcurve is assumed to be com- 
posed of statistically distributed shots of (possib ly) varying pro- 
files (cf. Lochner, Swank & Szymkowiak 1991 for detailed ref- 



erences and applications to Cyg X-l). Shot and distribution 
parameters are fit to various moments and statistics of the ob- 
served lightcurves. Time delays between hard and soft vari- 
ability are assumed to be due to differing shot profiles and/or 
shot distributi ons in different energy bands (cf. M iyamoto & 
Kitamoto 1989; N94; PF). Recently, PF have associated the en- 
ergy dependent shots with the 'energization' of a corona, which 
was parameterized by a series of equilibrium models where the 
energy input to the corona was a function of time. 

Does the data prefer one model over the other? The detailed 
structure of the PSD is likely to be more difficult to produce in 
the KHT model. Although the observed PSD can be modeled as 
a singly broken power law with an 'absorption' feature, it is dif- 
ficult to imagine a Compton corona configuration that would act 
as such a 'notch filter' . It is somewhat easier to imagine, how- 
ever, two or more preferred shot durations in the distribution of 
shot timescales. A suitable distribution of shot timescales could 
easily reproduce the type of PSD fits presented in §^|. 

Can the flux dependence of the PSD be reproduced? This de- 
pends upon whether the coronal radius decreases or increases as 
the observed flux decreases. The former possibility can be con- 
sistent with our 'sphere+disk' corona fits to the spectral data 



(paper I). The latter is more consistent with the ADAF models 
to the spectral data (paper I), as lower luminosity often implies 
a larger 'transition radius' to advective flow. The transition ra- 
dius in the ADAF models is the radius at which the flow transits 
from being geometrically thin, optically thick, and radiatively 
efficient to being geometrically thick, optically thin, and radia- 
tively inefficient. The larger this tra nsition radius, the smaller 
the overall accretion effici ency (cf. [Varayan 19%|; E sin, Mc- 
Clintock & Narayan 1997). Lower luminosities can be asso- 
ciated with lower overall efficiencies. The dependence of the 
coronal radius on observed X-ray lumi nosity is less const rained 
for our 'sphere+disk' coronal models ( pove et al. 1997 ; paper 
I). Larger coronae, however, can produce larger observed X-ray 
fluxes in these models. 

In the KHT model, the larger the corona the more attenu- 
ated the high frequency PSD will be. This will cause the PSD 
to apparently 'shift' towards lower frequency. Similarly, shot 
noise models usually associate the shot timescales with char- 
acteristic accretion flow timescales. What sets these timescales 
are usually not explicit in shot noise models; however, one ex- 
pects the timescales to increase for larger coronae, consistent 
with the expectations of the ADAF models. That is, we expect 
to see the lower characteristic frequencies associated with the 
outer accretion flow. Associating a smaller corona with lower 
luminosity, as can be fitted within the context of 'sphere+disk' 
corona models, would lead to trends opposite of the observa- 
tions. 

What is the expected relation between the PSD and the ob- 
served time lags for these models? If the coronal size decreases 
with decreasing luminosity, then the KHT model agrees with 
the observations. The shorter scattering path lengths of a small 
corona will produce characteristically shorter lags between the 
hard and soft photon variability. The KHT model also correctly 
reproduces the observed logarithmic energy dependence of the 
time lags. However, it is difficult for the KHT model to rec- 
oncile the lower frequency PSDs simultaneously with shorter 
time lags. Predictions for shot noise models are more depen- 
dent upon the specifics of each model. 

Energy dependent shots were first suggested by Miyamoto & 
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Kitamoto (1989) as the cause of the hard X-ray variability lags 
observed in Cyg X— 1. This model contained eight input shot 
profiles (four shot durations in two energy bands), and did not 
relate the timescales to specific disk physics. The magnitude 
and frequency dependence of the hard X-ray lags are adjustable 
by changing the duration of the hard shots relative to the soft 
shots. In describing observations of the 'very high state' tim- 
ing properties of GX 339—4, N94 related the shot durations to 
viscous and thermal instability timescales in an accretion disk. 
The shots considered in this work were also temperature de- 
pendent and became hotter as the shots progressed. This model 
reproduced the observed energy dependence of the lags, with 
the exception that it did not reproduce the 1.2-2.3 keV photon 
variability lagging behin d the 23-A.6 keV ph oton variability. 
As discussed in N94 and Mowak et al. (1997) , any shot model 
where the characteristic temperature or hardness of the shot 
smoothly increases as the shot rises will reproduce time lags 
qualitatively similar to those observed in Cyg X-l and GX 339- 
4. In the model of PF, the shot timescales are not related to any 
specific accretion timescales; however, the time-dependent shot 
hardness is related to a series of equilibrium Compton coronae 
models. There is a one-to-one correspondance between the shot 
amplitude and its hardness in this model. The PF model cor- 
rectly reproduces the magnitude and energy dependence of the 
observed time lags. 

We expect the following qualitative behavior for the mod- 
els of both N94 and PF. If we scale the shot duration to longer 
timescales, thereby shifting the PSD to lower frequencies, we 
also expect the timelags to increase, which is contrary to the 
observations. The only way to have a shot produce a lower 
frequency PSD and lead to shorter time lags is to alter the spec- 
trum of the shot as well as lengthen its duration. Whereas we 
do see changes in flux from GX 339-4, we do not see large 
changes in the best fit spectral parameters such as optical depth 
or coronal temperature. The observed energy spectra do not 
leave much room for varying the spectra of the presumed shots 
in order to fit both the high frequency PSD/long time lag data 
and the low frequency PSD/short time lag data that are ob- 
served. (It is possible, however, for the relative contributions 
to the spectra from a 'steady component' and a 'shot compo- 
nent' vary in such a way that the average spectrum remains rel- 
atively unchanged. The rms amplitude of the PSD significantly 
increases for the faintest observation, and therefore its average 
spectrum does have a greater contribution from the shots.) 

What is the expected coherence function for these various 
models? Both the KHT model and the PF model produce unity 
coherence, whereas the model of N94 produ ces coherence sub- 
stantially less than unity (Nowak et al. 1997). The KHT model 



value, a dip this low is very difficult to produce by summing 
different linear transfer functions with different phase delays, 
all the hard-to-soft variability phase differences for the individ- 
ual transfer functions are approximately the same magnitude as 
the observed phase lags. On the other hand, the nonlinear trans- 
fer functions of the N94 model are seen to lead to far stronger 
loss es of coherence, wit h typical coherence values being ps 0.1- 
0.3 (Nowak et al. 1997). Such coherence values are character- 



only considers static coronae, whereas the PF model uses the 
same exact linear Comptonization spectrum for each individual 
shot. We expect both of these situations to produce unity co- 
herence. The N94 model assumes both a sum of different shot 
spectra and each shot represe nts a nonlinear tran sfer function 
from soft to hard variability dNowak et al. 19971; Vaughan & 
Nowak 1997). 

Similar to Cyg X-l, GX 339^4 was seen to have near unity 
coherence over a broad range of frequencies, with a rolloff at 
w 3-10Hz. Unlike Cyg X-l ( [Nowak et al. 1998a| ), however, 
there was no evidence for a loss of coherence below w 0.02 Hz, 
but there was evidence for a coherence dip near 2 Hz (0.5 Hz for 
Observation 5). At first glance this seems consistent with both 
the KHT and PF models. We note, however, that the coherence 
function dips to as low as 0.9. Although still a large coherence 



istic of the coherence function observed at high frequency. (We 
previously have suggested that the observed loss of coherence 
above w 3- 1 Hz may be related to nonlinear effects due to flar- 
ing activity feeding a corona on dynamical timescales; N owak 
etal. 1998a.) 

The question put forth in |Vaughan & Nowak (1997) was: 
why are observed coherences so close to unity? With the ob- 
servations presented here, and in light of the discussion of §^ 
and equations (Q) and (Q), we might wish to modify this ques- 
tion to: when the coherence function is so close to unity, why 
is it as low as observed? We argued in §|] that if there truly 
is a relationship between the time lags and the coherence, then 
there must be a mix of processes with a broad range of time 
lags. The fits to the phase lag and coherence data suggested one 
component with near zero lag between hard and soft variability. 
The fits further suggested that the higher frequency Lorentzian 
component exhibited characteristically longer time lags (at a 
fixed Fourier frequency) than the lower frequency Lorentzian 
component. This possibility was not considered in the models 
of N94, KHT, or PF. However, as the KHT model postulates 
only one 'transfer function' from soft variability to hard vari- 
ability (i.e., diffusion through the corona), it is more difficult 
to reconcile this model with the conjecture that phase lags and 
coherence are created by a sum of independent processes. Shot 
noise models, however, can introduce multiple, albeit currently 
unknown, independent processes in the manner suggested by 
the fits to the data. 

6. SUMMARY 

We have presented timing analysis of a series of 8 RXTE 
observations of the black hole candidate GX 339-4. On long 
timescales, there is evidence of a 240 day periodicity in the 
ASM lightcurve. This is not a strictly periodic feature, but 
is probably more like a 'characteristic timescale'. Such a 
t imescale is consistent with warping and precession timescales 
(Pringle 1996; Maloney, Begelman & Pringle 1996); however, 
as we discuss in paper I, the evidence points more towards coro- 
nal size variations than to inclination effects. 

In terms of characterizing the short timescale variability, we 
see that the observations break up into two sets. The seven 
brightest observations span a range of two in observed 3-9 keV 
flux and all have comparable timing properties. These observa- 
tions all show evidence of an ps 0.3 Hz QPO. The properties of 
this QPO are not strictly steady; however, there are no obvious 
correlations with the flux of the source. 

The faintest observation, which is 5 times fainter than the 
brightest observation, had: larger amplitude variability, charac- 
teristic PSD frequencies that were a factor of three lower than 
observed for the other PSDs, and showed shorter time lags be- 
tween hard and soft variability. This latter feature, albeit with 
a great deal of scatter about the trend, was also mirrored in the 
flux dependence of the time lags observed for the brightest ob- 
servations. 

All PSDs were reasonably well-fit by the sum of three fit 
components: a power law, a low frequency Lorentzian, and a 
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high frequency Lorentzian. We further suggested that the ob- 
served Fourier frequency-dependent phase lags and coherence 
could be explained by summing these three components with 
differing intrinsic time delays between hard and soft photon 
variability. 

We discussed all these possibilities in light of various theoret- 
ical models. The short time delays for the lowest flux observa- 
tion ap pear to be more in ag reement with 'propagation models' 
(KHT; |^Iowak et al. 1998b ) if the coronal size decreases with 



decreasing flux. If, however, as suggested by ADAF models the 
coronal (i.e., advection dominated) region grows with decreas- 
ing flux, the trends observed in the time delays are counter to 
the theory. Conversely, the low flux/low frequency PSD ob- 
servations are more in agreement with the models where the 
coronal region grows with decreasing luminosity. Shot noise 
models are more likely able to be adapted to explain simulta- 
neously the phase lags and coherence as a sum of independent 



linear components. 

None of the models currently address the complicated 
nonlinear processes that may be occurring on dynamical 
timescales. Such processes, which could be the 'flares' that 
energize the corona a nd/or the observed r adio-emitting outflow 
( |Fender et al. 19971; [Corbel et al. 19971; H annikainnen et al. 
1998; paper I), could be the cause of the strong loss of coher- 
ence seen at high fre quency in both GX 339^1 and Cyg X-l 
( |Nowak et al. 1998tj ), as well as the cause of the 'flattening' 
with increasing photon en ergy observed at the high-frequency 
end of the Cyg X-l PSD ([Nowak et al. 1998a[). 
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